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A non-perturbative scheme, based on the functional generalization of the Callan-Symanzik equa- 
tion is developed to treat the Coulomb interaction in an electron gas. The one-particle irreducible 
vertex functions are shown to satisfy an evolution equation whose initial condition is given by means 
of the classical action and the final point corresponds to the physical system. This equation is trun- 
cated by expanding it in momenta and excitation energies, leaving the electric charge as an arbitrary, 
not necesseraly small parameter. Exact coupled partial differential equations up to first order in the 
frequencies and excitation energies are derived. The numerical integration of these eqations is left to 
a later stage. Nevertheless, in order to demonstrate the relation with the perturbation expansion the 
one-loop Lindhard function and screening are reproduced in the independent mode approximation 
of the evolution equation. 



I. INTRODUCTION 



In the many electron problem the Coulomb interaction is usually taken into account by a partial resummation of the 
perturbation expansion. A possible non-perturbative approach, based on the concepts of the renormalization group 
and the Callan-Symanzik equation is proposed in this paper. Originally viewed as a physically motivated partial 
resummation of the perturbative contributions [0, the renormalization group strategy has lately developed into a 
general method to treat strongly coupled systems ||- ||. The advantage of this method comes from the functional 
treatment of the renormalization group flow equation where the only approximation is the truncation of the effective 
action according to the gradient expansion. This formalism, originating from the Wegner-Houghton equation [0, 
allows us to follow a large number of coupling constants and makes the approximation rather flexible. Although the 
Callan-Symanzik || equation is not a renormalization group equation because it traces the changes of the dynamics 
as the mass is varied, its appearance and universality, assumed that the mass is the only scale in the UV regime, show 
similarity with the renormalization group method. 

The traditional application of the renormalization group method for fermions with finite density JTT[ represents 
a promising approach for the treatment of the Coulomb interaction in an electron gas. We present in this paper a 
version of the functional renormalization group method for systems with Fermi surface which allows for the systematic 
improvement of the scheme followed in these works. The complication of the blocking in fermionic systems comes 
from the fact that the successively eliminated modes should be concentrated more and more at the Fermi surface as 
we follow the renormalization group flow. In the zero density case the infrared limit consists of a single homogeneous 
mode. But the modes eliminated by the usual blocking in the infrared regime of fermionic systems with finite density 
populate the Fermi surface. This poses a problem for the implementation of the usual blocking procedure where the 
interactions are represented by fields. On the one hand, the momentum exchange between electrons is kept finite, of 
order of in the infrared limit, which requires the evolution of these coupling strength. On the other hand, the scalar 
fields are reduced to the point p = in the infrared limit and their coupling strength at the scale hp is already frozen 
out. In other words, the momentum modes of the field variables, responsible of the Coulomb interaction can not be 
eliminated in a sequential step-by-step manner during the blocking and one can not avoid non-local interactions. 

In order to overcome this difficulty we do not order the modes for their elimination according to their energy. 
Instead, our blocking will eliminate the quantum fluctuations in the increasing order of their amplitude In this 
manner we can keep the fundamental requirement of the blocking, namely that it starts with the more weakly coupled 
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modes and treats the more strongly coupled ones later in the renormalization group flow, without paying a high price 
for the presence of a Fermi surface. 

The strategy of this generalized renormalization group method can briefly be summarized as follows. First, one 
introduces a parameter M in the Hamiltonian in such a way that H — > Hm where Hm is almost diagonal in the free 
particle basis for large M, i.e. the interactions are suppressed as M — > oo. The next step is the computation of the 
evolution equation in M for the effective action, the generator functional of the one-particle irreducible (1PI) vertices. 
The renormalization group improved approach consists of a systematic expansion in the vicinity of the current system, 
namely considering Hm as the non-perturbed Hamiltonian and using AHm = Hm-am — Hm as perturbation to 
obtain the changes in the effective action as M is decreased by AM. Notice that the exact evolution equation can 
be obtained already in the leading order of the perturbation expansion as AM — > because AHm = O(AM), in 
other words, AM/M is a new small parameter to obtain a differential equation for the M-dependence. Once the 
evolution equation is known, its integration from the perturbative initial condition of the system from large M down 
to M = produces the desired 1PI vertices. The evolution equation, being a differential equation for a functional, 
becomes useful only after projecting it onto a restricted subspace of functionals. In this subspace the effective action is 
characterized by ordinary functions and the evolution equation is converted into a set of coupled differential equations 
for them. 

This method is based on the following two general assumptions: (i) The dependence of the 1PI functions on the 
control parameter is supposed to be differentiable. Notice that this condition, the differentiable "turning on" the 
Coulomb interaction is weaker than the assumption needed for the applicability of the perturbation expansion, the 
analyticity for arbitrary values of the control parameter, (ii) The possibility of using a simple ansatz for the effective 
action to incorporate the physics of the model. This step turns the evolution equation, a functional differential 
equation into a set of coupled differential equations. With these assumptions the perturbation expansion around the 
Gaussian system is avoided, the 1PI functions are obtained by solving differential equations rather than computing 
Feynman graphs. 

The parameter M controls the amplitude of the quantum fluctuations. In a relativistic system the mass plays such 
a role and the simplest realization of the blocking in the amplitude is based on the modification m 2 — » m 2 + M 2 of 
the Hamiltonian. The resulting differential equation is the functional generalization of the Callan-Symanzik equations 

i- 

In this paper we present the first step only in the implementation of these ideas in the framework of the Coulomb 
gas. We consider 1PI functions with the following restrictions: (i) The dependence on the external energies and 
momenta is kept up to 0(ui) and 0(£) for electrons where 10 is the frequency and £ is the energy counted from the 
Fermi surface and up to O(oj 2 ) and 0(p 2 ) for photons, (ii) There are no more than two external electron legs. The 
number of photons is not restricted. 

We arrive at a set of coupled differential equations in this manner for the functions parametrising the 1PI functions. 
The numerical integration in the control parameter yields the solution of the Coulomb problem for the assumed ansatz 
for the 1PI functions. There is no problem in principle to generalize the scheme and include 1PI functions with more 
external electron lines and more complicated energy and momentum dependence. 

In order to demonstrate the relation between the perturbation expansion and our scheme we integrate out the 
evolution equation in the independent mode approximation. In this approximation the system remains unchanged 
during the evolution, meaning that the modes are taken into account independently from each other. This corresponds 
to the one-loop approximation of the perturbation expansion. The one-loop Thomas-Fermi screening and Lindhard 
function are reproduced in our scheme without considering Feynman graphs. 

We derive the evolution equation for the electron gas with Coulomb interaction in section II. The resulting functional 
differential equation is projected on the leading order gradient expansion ansatz (expansion in u>, £ and p) for the 
effective action in section III. In order to check the consistency of this scheme with the standard perturbation expansion 
we solve the evolution equations analytically in the independent mode approximation in section IV. Section V contains 
the conclusions. Some technical details of the functional derivation and the gradient expansion are shown in the 
Appendices. 

II. EVOLUTION EQUATION 

We shall study spinless electrons interacting by means of the Coulomb potential and described by the Hamiltonian 

r r h 2 1 r e 2 

H= d 3 x^(x) A-p tp(x)+ d 3 xd 3 y^(x)^(y) — ; ^(yU(x). (1) 

J [ 2m J 4tt\x - y\ 
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where the field ip corresponds to electrons and \x is the chemical potential. In the path integral formalism we introduce 
an auxiliary field u(x,t), taken formally as the temporal component of the photon field, write the action in terms of 
the Grassmannian ip, ipi anc j rea i u as 



t „i - 



ihdt H A + jj, + eu x 

2m 



in the functional integral 



V[i^]V[i)]V[u]e 



(2) 



(3) 



The x, t or p, oj arguments of the field variables will be denoted by the subscripts x and p, respectively. The symbols 
x and p stand for the spatial components only if written otherwise. In a similar manner, we shall use the notation 
f x = fdtd 3 x = VT, L = J divd 3 p/(2n) 4 , and f x = f e iq - x - iut f q . The field u x can be considered either as an 
auxiliary field, responsible for the Coulomb interaction or the temporal component of the photon field in the Feynman 
gauge after both the transverse and longitudinal components of the spatial photon field have been suppressed. The 
homogeneous mode of the photon field, 



u 



1 

VT 



(4) 



plays a special role because of the absence of 0( u o) terms in the action. In fact, the homogeneous component of 
the photon field uq contributes to the chemical potential modifying the electron density and the inhomogeneous field 
components induce polarization effects only. Another way to see this is to note that the integration over the uq 
configuration induces the constraint f ip x ip x = and wipes out asymptotic states with non-vanishing electric charge. 



Actually, the homogeneous photon field configurations in QED, A fi (x,t) 



are not to be integrated over in the 



path integral of a non-confining gauge theory ]12| . Thus the integral measure T>[u] does not involve the integration 
over the homogeneous mode i.e. uq = 0. 

We add a quadratic term to the action, S — > Sm = S + S s , 



S.[i>,^,u]= / /(M 2 )V4 U-A + mW 



2m 



M 



(5) 



in order to control the quantum fluctuations. The function f(M 2 ) is chosen in such a manner that /(0) = and 
liniM— »oo f(M 2 ) — oo, in order to suppress the fluctuations of the fields when M — > oo. In practice we shall use the 
simple form /(M 2 ) = M 2 /m 2 . Since we are interested in a system of particles at finite density, the chemical potential 
fi is non- vanishing. This particular form, (||), was chosen to control the strength of interactions by changing the 
weight of the quadratic, non- interacting part of the action. The gradual decrease of S s "turns on" the interactions 
and the integration of the evolution equation in the control parameter allows us to solve the interactive system. 
The suppression action can be chosen to be gauge invariant, but we are satisfied with this simpler form because 
the restriction for the Coulomb interaction only, the suppression of the spatial component of the photon field alone 
removes the gauge symmetry. 

The action Sm determines the dependence of the microscopic physics on the control parameter M. In order to trace 
down the role played by M we consider the generator functional of the 1PI vertices, the effective potential defined by 



where / • g = J x f x g x , 



(6) 



(7) 



and 



6W M \j,jKJ] 







sw M [j,j^ 


J] 








J] 



SJ X 



( U x)j,jf,J = U x- 



(8) 
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We use the same notation for the quantum field and its expectation value. 

We first derive a functional differential equation in M 2 for the effective action, 



(IT 



M 



d{M 2 



(9) 



whose integral, 



To — fjif — 



M- 



d(M') 2 F[r M ,;M'] 



(10) 



yields the effective action of the model we want to solve in terms of the initial condition Tm[iP, i/j t, u] = Sm[iP, > u ] 
imposed when the M 2 exceeds every energy scale of the model. 
The evolution equation of the effective action can be written |§ as 



d M 2r M [ip,^,u] = d M *W\j,j\ J], 



(11) 



where 8m 2 = d/dM 2 according to (^|) and (||). Since we know the control parameter dependence of the bare action 
we have immediately 



W M 



(12) 



where /' = 9m 2 /, £ — —h 2 A/2m — /i measures the kinetic energy operator measured from the Fermi level and Tr 
denotes the functional trace, TrA — J A x>x . 
We introduce now the notation 



(1) _ 6_F[f] (2) _ S 2 F[f,g] 
/ ~ Sf ' f9 



Sf6g 



(13) 



to relate the second functional derivatives of W and f 



/ Sip z 6u z \ 

Sipi Sipi Sipi 

Sipl 8ip z Su z 

i^x &ipx 

\ Sipl 5ip z Su z I 

\ 8u x 8u x Su x ' 



/ f (2) f (2) f (2) \ / ^(S 

ViVy ViVy Vi^t 



f(2) f(2) p(2) 

i>xll>y lp x 4>l 



f(2) 

\ «xV>y 

The evolution equation (Jll|) goes over to 



f(2) f(2) 

«x-0j 



d M *f M \^^,u] = -ihf'(M 2 )Tr[EWj*j t ] - f \M 2 )Tt[£^ >] + yTr^j 2 ]] - ^Tr[ U 2 ] 



-W (2) 


TyTz 


JyJ z 


w {2) 

3y3z 


-w {2 \ 

jyjz 




wf) 


J yJz 





if i. 



= 5 x .z ■ 1 



(2) 



(14) 



(15) 



In order to remove the bilinear terms in the fields on the right hand side we write out explicitly the artificial suppression 
terms in the effective action, 



f m [V' t » "0, u] = Im u] 



f(M 2 )^ x £^ x + ~-M 2 u 



1 



2„,2 
x ] ! 



(16) 



which leads to 



d M 2Y M [^^,u] = -ihf'(M 2 )Tr 



E^-M 2 



-Tr 



a closed evolution equation for the effective action, with 

/-£ \ 
M 2 = £ . 
V M 2 J 

Two remarks are in order about details not shown explicitly in Eq. ([l7|): 



J, J 



(17) 



(18) 
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• One should bear in mind that the right hand side is UV divergent as it stands because the suppression of 
the amplitude of the fluctuations is not a regulator. In order to have well defined expressions we introduce a 
cut-off A for the momentum integrals. This regulator will be kept finite because our non-relativistic model is 
non-renormalizable . 

• We are working in real time, and the quadratic part of the action must be supplied with an infinitesimal 
imaginary part. We shall use the electron propagator 

G(p, u) = P- = { eik /- p) + e( T fef) ) eW , (19 ) 

where 77' — > + . Since the chemical potential dependence is not always continuous we have to make sure that Uf 
does not evolve with the control parameter M 2 in order to avoid the appearance of singularities in the evolution 
equation. We shall achieve the M-independcncc of the Fermi-momentum by the choice S p — + sign(p — ftp). 
In the case when the chemical potential dependence is continuous then schemes with an M-dependent Fermi- 
momentum are possible but their evolution equation appears with an additional term, d^i/dM times the density 
in©. 

III. GRADIENT EXPANSION 
A. Effective action 

In order to convert the functional differential equation ( |l7| ) into a set of coupled differential equation for the running 
coupling constants which are the parameters of the effective action we introduce an ansatz for r^f ip, u]. We first 
separate the fluctuations of the photon field from the homogeneous background, 

u x = Uq + T]x, (20) 

where rj p= Q = 0. Note that once an ansatz is introduced for the effective action we may explore its structure for uq 
even if homogeneous photon field configurations are not allowed in the path integral. Thus uq ^ will be used in the 
effective action whose form will be written as 



I 1 

iX (u x )d t - Xi(ux)£ + U x {u x ) + id £ + -Y (u x )[d t u x y + ~Yi(u*)(V 'u x f 

^Z Q (u x ){d t u x ) 2 - iz 1 (u :E )(Vu :I ,) 2 - U 2 (u x ) +ieu 2 x 



ipx 

(21) 



where 5g — + sign£ and the Xj, Yi, Z% and Ui are polynomials in the field u x . This funtional form is given in the 
framework of the gradient expansion which is an expansion in the space-time derivatives of the photon field truncated 
at 0(d 2 ). We write the effective action in terms of the 1PI functions in momentum space, 



J pi,P2 •>pi,pi,q 

+ \l ^Wn + l f r£iZu^k^ u ^ + --- ( 22 ) 

•'91,92 ^P1,P2,91,92 

where the superscripts show the type of the external legs and the subscripts give the corresponding momenta. One 
finds by using this form that the gradient expansion corresponds to the Taylor expansion of the 1PI functions in 
the photon energy-momentum. In the present work the electronic 1PI functions will be considered by expanding in 
the frequency and the kinetic energy, dt and £, respectively. The dependence on the Grassmannian fcrmionic field 
variables is kept quadratic in the present approximation which reduces the argument of the local coefficient functions 
X, Y, Z and U to the photon field u. 

The potentials U\(u) and U2(u) are space-time independent, indicating that they do shift the contour of the 
frequency integrals. The finite lifetime effects might come from the imaginary part of the electron self energy. This 
appears in 0(uj 2 ) and will be ignored in the present work which includes 0(lu) only. 

The infrared divergences of the perturbation expansion generate singular behaviour at w = p = 0. This fact 
introduces some complications. In particular, the photon 1PI functions are not continuous at this point and their 
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limiting value depends on the way this point is approached, as the remnant of gauge invariance, the transversality 
of the photon self energy. This feature introduces an important difference between the 0(rf) and the 0(rj 2 ) level 
solution of the evolution equations. We have already mentioned after Eq. ^ that the dynamics of the homogeneous 
mode of the photon field is different from the non-homogeneous ones. To understand better the way this appears in 
our scheme one must remember that the solution of the complete functional evolution equation for the effective action 
generates the exact 1PI functions, cf. Eq. (^). As we shall show below, the truncated gradient expansion ansatz for 
the effective action reduces the evolution equation into a hierarchy of approximate relations between the 1PI functions 
because only a polynomial of finite order in the external momenta is retained for the 1PI functions in the otherwise 
exact equation. The integration of the contributions of order 0(i] ) to the evolution equation reproduces the value 
of the 1PI functions at u> = p = in this approximation. For example, the integration of the terms 0(r] 2 ) yield the 
photonic two point 1PI function at arbitrary energy-momenta T^"" (u%,Pi, u>2,p2), but 

iimrr( w , P ,- w ,- J) )^r;r(o 1 o,o,o). (23) 

We take such a non-continuous energy-momentum dependence of the effective action into account by replacing the 
potential Ui{Ux) in by two different functions, 



2[x) \VTU 2 (uo) 77 = 0. 1 ' 



The first line on the right hand side, representing the left hand side in Eq. (|23|), will be evaluated by taking the limit 
u -> before p -> 0, cf. Eq. @. 



For large enough M, when the quantum fluctuations are negligible, we have the initial condition 

X = ft, X 1 =Z 1 = 1, Y = Y 1 = U 2 = U 2 = 0, Ui = eu, (25) 
up to 0(M~ 2 ) corrections because Tm[^ , ip, u] — S[^, ip, u], see after Eq. (|l0|). 



B. Expansion in the photon held fluctuations 



In order to determine the evolution equations for the coefficient functions we write the photon field as in Eq. ( pOj ) 
and expand Eq. (|l7|) in powers of r] x . The left hand side is 



T M [i/>1, 4>, u] = T M [^, i>, u Q ] + J r) z T$ Uz [^,i>, u Q ] + ^J ^ VzVz'T®^, uo] + 0( V 3 ), (26) 

when r) ^ 0. The zeroth order term in r\ gives 

r M [^, V, u ] = f il>l q %A q - VTU 2 , (27) 

Jq 

where Uq is an arbitrary parameter, introduced in order to explore our functional differential equation (JTt]) , 



A„ = X u; Q - X\Ea + J7i, 



(28) 



the coefficient functions are taken at uq, and £ q = q 2 /2m — \i. Note that even though we have u — uq in Eq. ( |27| ) 
we are dealing with the r\ ^ part of the effective action because the functional Taylor expansion coefficients are 
obtained in the limit rj — > 0. The contributions to the first and second order in 77 can be written as 



•'91,92,93 



and 



L 



Vqi Vq 2 ^qi ^ ~ Y UJ qi U q2 - Y X q\ ' ^2 



91,92,93,94 



93 "91+92+93,0 



3 9l+92+93+94,0 



VqV-q 



Y io 2 + Y iq 2 + C/ 2 (2) 



(29) 



(30) 



by means of the functional derivatives ([All) and (A3) recorded in the Appendix. 



G 



On the right hand side of the evolution equation ([D]) we find the second functional derivative matrix f jj^ = 
Td + r n d + Tfl, where Td Pl , P2 is vanishing for pi + P2 ^ 0, T n d = 0(rp) and Tn comprises the 0(rj) and 0{rj 2 ) 
contributions. The O(rj ), diagonal piece is obviously 





( ~ G V2 













G P1 







V o 








3 Pl+P2,0i 



(31) 



where 

G; 1 = A p -f(M 2 )S p + iS £p , 

Dp 1 = Z Q uj 2 p - Z lP 2 - C/f 3 - M 2 + ie. (32) 

Note that the U\ of Eq. ( g^ ) does not influence the imaginary part of G and the Fermi level remains independent of 
the photon field. This is in agreement with our introduction of uq in the effective action as a formal parameter to 
help to uncover the content of the functional equation (|l2| ) rather than a dynamical variable. 
The remaining 0{rp) contributions are easily found, 

/ o y,^ u * 
r„ d = £^ u | (33) 



the matrix elements being listed in Eq. ( A5). The terms 0(rj) and 0(i] 2 ), parametrized as 

1/1+1/1 (T V> t '»/' t rr^ u 



a W o-V-V-t a 4>u ] +0{t 1 3 ), (34) 



are explicitly presented in Eq. (A9). 

The inversion of r*- 2 -* is carried out by expanding in the non-diagonal pieces, 

W< 2 > = W n (2) + V + V W^ 2 > + V + V +■■■ (35) 

' / y "oi 1 / , "ni,a 2 i / j " 01,02,03 ~ / y "01,02,03,04 1 V.""/ 

Oi Qi,0 2 01,02,03 01,02,03,04 

(2) —1 

where the summation is over the indices ctj — nd or fl and Wq = T d denotes the diagonal part of the propagator. 
The non-diagonal contributions are collected in 

"qi vt "i ' 

<!o2,o 3 = ^o (2) roi^ (2) r Q2 ^ (2) r Q3 iy (2) , 

^i 2 L,03.o 4 - -W^T ai wPT a2 W^T a3 W^T a M 2) - (36) 



In the last step we substitute the gradient expanded form, ( |21| ) , and identify in (jTTj) the contributions of the different 
powers in the electron field and the fluctuations in order to get the evolution equations for the coefficient functions. 



C. Homogeneous photon field 



The zeroth order in the fluctuations is obtained by setting 77 = on both sides of (17). On the right hand side, in 
particular in Eq. (p5[) this amounts to neglect terms containing Tr. The rest involves the expressions 



Tr 



(lf)-M 2 ) 



3,3' 



= -T:[£WW] = 



£pG p 



P J P < 



c v Kj v lj - P , P ,u v 



(37) 



and 
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Tr 



./../ 



A 2 )] 



= / [D p - D^l p ] + [ D 



(38) 



P ~P f ,P 



The left hand side of the evolution equation (|l7| ) includes the expression (|7|) and the use of ( j32| ) and ( A5) gives the 
evolution equations 



9m 2 ^2 = ihf 



— ih 



p x o^ P - (Xi + f)£p + Ui J p z 0bJ 2 _ Zip 2 _ jj_ 



r(2) 



M 2 



(39) 



&m 2 [XaU q — Xi£ q + Ui] 



itl 



A q 2) + Y u; 2 + Y lP 2 



2 J P (Z u 2 - Z lP 2 - U^ 2) - M 2 ) 2 



A 



(i) 



-A* 1 , 

2 q Jp (Z u 2 -Z^ + q^-uP -M^(X oj p -(X 1 +f)£ p + U 1 ) 



-ihf'AW 



^p-™-—p 



p (x lu p - {x x + f)e p + Ux) 2 (z^ 2 p+q - z x (jp + q y - uf ] - m 2 ) 



(40) 



We introduced the symbol to indicate that the right hand side is supposed to be written by expanding in the 
variables lu, £ and p and by retaining only the orders shown on the left hand side. 

As we shall see below, the discontinuity discussed in Eq. ( |23| ) affects only Eq. (|39|). Eq. (^) gives the evolution of 
the coefficents Xq, Xi and U%. In order to solve these coupled differential equations one needs other equations for the 
Zi% YSs and U 2 . 



D. Fluctuations 



In order to find the evolution of the coefficient functions we substitute the gradient expanded form, (pT|), and 
identify on the right hand side of Eq. ( |l7|) the contributions of different powers in ipty and ij. Some expressions of 
the contributions of the individual terms on the right hand side of Eq. (^) are shown in Appendix |b| 

It is easy to isolate the terms O((t(j^ip) ri 2 ), 



~ (d M 2Z Q u, 

ihf 



d M iZiq 2 



f A {1) A {1) 



p (X u;p ~ {X x + f)£ P + Ux) 2 (X uj p+q - (Xx + f)£ p+q + Ux) 



+ i 



£ P {A 



(2) 



Y LU 2 + Yxq 2 



ih 



„ (X u p -(Xx + f)£ p + Ux¥ 



ih 

T 



gg (^ + ^) + zr(g 2 +p 2 ) + ^ 
(Z cj 2 - Z lP 2 - C/ 2 (2) - A/ 2 ) 



(2), 



[Z^> (co 2 p + J\ + LUpLOg + Z\ L > (p 2 + q 2 +p-q) + U, 



r( 3 )l2 



P (Z oW 2 - Z lP 2 - tff - Af 2 ) 2 (Z c 2 +9 - Zi(p + </) 2 - U?> - M 2 ) 



(2) 



(41) 



The discontinuity ( |23] ) appears here as the non-commutivity of the integration and the limit u> q , q — > in this equation. 
The singularity at g = is caused by the fact that the order of the singularity in ui p is 2 or 3 when q — or q ^ 0. 
Since we discuss the 0(rf) contributions we have to keep q ^ 0. The g-independent part of the right hand side, the 
extrapolation from q ^ to q = will contribute to rather than according to Eq. (p4|) . 



The terms O^tprf) on the left hand side of Eq. tUTT) is 



J M 2 



9l-?2 ^92 



91,92 



A, 



^92 



(1) 



(42) 
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Since we have q± ^ there will be no singular coincidence of the poles on the right hand side of the evolution equation. 
As a result the terms 0(w° 2 £g 2 ), C(w 92 f° 2 ) and G(u)® 2 £ q2 ) on the right hand side of the evolution equation of order 
0{^ipr]) reproduce the first derivative of the right hand side of Eq. ( ff(i| ) with respect to wo- 

In the order 0(ip^ipr] 2 ) we find 

f + " (2) (2) (2) 

^<9 M 2 / Vq^T-q! -q 2 - q3 i>q 3 ^0 W 93 ~ ^1 £q 3 + U \ ~ Y QU qi UJ q2 - Y iqi ■ q 2 (43) 

on the left hand side of the evolution equation. The terms 0{ui q3 ), 0(£ q3 ) and 0(1) give the second derivatives of 
Eq. ([io| ) and can be omitted on both sides of the equation. The remaining contribution in 0{u> qi uj q2 ) and 0(q\ -92), 
obtained for (73 = is 



/ VqiVq2^- qi - q2 - q3 ^q3 l>Q W «l W «2 + Y l1l ' 92] 

•'(71, <72,03 



'01,92,93 

S th (Tr [W fljJ »] -Ti-Tz- Ti) - f(M 2 ) (Tr [Wfl^j] + B 2 + £ 3 + B 4 ) , (44) 

where the quantities on the right hand side are given in Appendix [b| When the coefficients of the field variables with 
the same arguments are identified on both side we find two differential equations, one for Y$ and another for Y\. 

The form of these and the other differential equations for the coefficient functions can simply be understood in the 
following manner. We start with the one-loop expressions for the 1PI vertex functions using the non-Gaussian parts 
of the effective action T as interaction. The infinitesimal increase of the parameter M 2 — > M 2 + 5M 2 modifies the 
electron and photon propagators as 

G -> G - SM 2 f'{M 2 )G 2 , D -> D - SM 2 D 2 , (45) 

respectively. The change induced by this transformation in the one-loop expressions for the 1PI vertex functions is the 
right hand side of our evolution equation up to the factor 5M 2 . At this stage we recover the perturbation expansion 
based on the Feynman rules. The usual accumulation effect of the renormalization group strategy occurs when the 
differential equations are integrated out. By making n steps M 2 — > M 2 + SM 2 the integration sums up the n-loop 
contributions to the corresponding 1PI vertex function. 

In this manner the integration of the evolution equations Eqs. (^), (f4(i|), ( |4l| ) and ( |i4| ) with the initial condition 
( p5| ) provides the effective action, the 1PI functions in the given approximation of their dependence on the external 
momenta. The successive inclusion of terms O(('0^'0)' i ) or higher order derivatives in the ansatz ( ^l| ) for the effective 
action yields more complete solutions. 



IV. INDEPENDENT MODE APPROXIMATION 



The renormalization group strategy is based on the accumulation of the effects of the modes which are eliminated 
sequentially. This corresponds to the integration of the evolution equations where the right hand side evolves with 
M 2 . We can easily obtain an analytic approximation for the evolution by ignoring the M 2 -dependence on the right 
hand side of the evolution equation. This amounts to the independent mode approximation which reproduces the 
one-loop 1PI functions in the given order of the external momenta. 



A. Homogeneous photon field: Hartree-Fock approximation 



As already mentioned in paragraph III A, the 1PI functions are non-analytic at u> = p = 0. We shall first show 
that this implies that U2 does not evolve. The right hand side of Eq. (|l7| ) is obtained by using an M 2 -independent 
effective action (|2l|), given by Eq. (Ba). As far as the 77 = case is concerned we find 



d M *U2 = -ifif 



l p hu p - (1 + f)£ p + eu 
After the frequency integration of the first term, 



ih 



I p p 2 + 



M 2 



lim 



d 3 p 



f 



du 



ri^oj (2tt) 3 p J 2n hujp - (1 + f)£ p + eu + iS p 



^ £piQ(—£ p ) = —i 



(27T) 



307r 2 m ' 



(46) 



(47) 







the right hand side turns out to be uo-independent. Since we are interested in the field dependent part of the effective 
action only, Eq. ( |46|) gives Om 2 ^ — 0. Polarization effects do not appear for a homogeneous field. 

The non-analyticity of the effective action at uo — p = is due to the photon sector only. The electron sector can 
be studied in the homogeneous photon field approximation, rj = in Eq. (|2^). The identification of the contributions 
0(V> f iO in Eq. @, the use of Eqs (§f) and @ taken at the one-loop approximation and the initial condition 
Eq. (pit ) lead to 



M 2 [X io q - X x £ q + Ui] E zhf'e 2 f 

J V 



f 



(hcop - (1 + f)£ p + eu Q ) 2 ((p + q) 2 + M 2 ) 



+ihe 2 I =■. (48) 

J P (hcjp - (1 + f)£p + eu ) ((p + q) 2 + M 2 f 

The frequency integral in the first term appearing in Eq. ( [48| ) vanishes. The second term gives 

_ i d? P e(-Sp) ^ p 2 d P 

J (2tt)3 (( P + q ) 2 + M 2 ) 2 2tt 2 J ( P 2 + q 2 + M 2 y~W 1 ' 

In order to identify the O(oj q ), 0(£ q ) and 0(1) terms with those appearing in the left hand side of Eq. (ff8l), the last 
expression is expanded in powers of £ q , 

- 2 < k ' p2dp f i- ,^:*i +M :L s q+ o(s?)). , 5 o) 



2^ 2 Jo (p 2 + k 2 F + M 2 ) 2 - 4p 2 k 2 F \ {p 2 + k 2 F + M 2 ) 2 - 4p 2 k 2 F q K q 
The identification of the different terms in Eq. ([48]) leads to the evolution equations 8m 2 Xq = 0, and 



d M 2 Xi = 



-2me 2 f kF p 2 dp(-p 2 + k 2 F + M 2 ) 



{{p 2 + k 2 F + M 2 ) 2 - Ap 2 k 2 F f 



„ rr ^ r «U 

dM2Ul " 2^io (p 2 + k 2 F + M 2 ) 2 -4p 2 k 2 F [bL) 

The right hand side of these equations is ito-independent, giving X 0M 2 = h, Xi m 2 ( u o) = ^ijf and U\^[uq) = 
Ui,m 2 = eu a = with 

me 2 f 1 , M 2 1 1 



d M iX 1M2 - -— |^3- In 4jfc 2 +M2 + M4A .2 +M2) + ( 52 ) 

In order to obtain the expression for X\ m 2 i one nas to integrate this diffential equation over M 2 from M 2 -> oo to 
e — > 0. This yields 

Trie 2 \ f € \ 

X ^ = ~^T F [ 2 + l^W F J >• (53) 

This coupling constant represents the ratio of the bare to the effective mass of the electron. It diverges as e — > 
meaning that the electron effective mass vanishes. This is a well known result from perturbation theory ]l5fi and 
corresponds to the Hartree-Fock approximation. Note that the divergence in the perturbation expansion is also 
logarithmic and appears as lhrifc^fc F ln(fc — kp)/2kF- The singularity appearing when k — > k F corresponds in the 
M-parameter space to M — * 0. 



B. Fluctuations: Thomas-Fermi mass 



The reason to go beyond the homogeneous photon field approximation is the need of the kinetic energy for the 
photon to reproduce the desired non-analycity in the photon self energy at u> = p = 0. The terms 0(rj 2 ) in the 
evolution equation determine the inverse photon propagator, 



l - (d M 2Z LO 2 q - d M 2Z lQ 2 - d M 2U {2) ) S 



d 3 p f dujp £ p 



ihf'e 2 / t- I —L zE (54) 

(27r) 3 J 2ir (hujp - (1 + f)£ p + eu ) 2 (hujp +q - (1 + f)£ p+q + eu ) 
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The frequency integral yields 



-Q(-s p ) + e(-^ 



p+q) 



which can be written as 



h (htu q - (1 + f)(£ p+q - £ p ) + i(S p+q - S p )f 
Q(-£ p )Q{£ p+q ) Q{£ p )Q{-£ p+q ) 



K P \ {hw q - (1 + f){£ p+q - Sp) + i5f (hu q - (1 + f)(£ p+q - £ p ) - iS) 2 $ ' 
remembering that 5 P = sign(£ p )S. Performing the change of variable — p — q — > p in the second term, one finds 

\(d M ,Z Q uj\ - d M2 Z iq 2 - d M ,U (2) ) 55 f'e 2 J j0-Q(-£ p )[l - Q(-£ p 



(55) 



(56) 



■'P+q) 



■■p+q 



{hu q - (1 + f)(£ p+q - £ p ) + iS) 2 {hu q + (1 + f){£ p+q - £ p ) - i8) 2 



(57) 



The dependence on u q drops out as q — » 0, giving Z$ = Ti. 

In the following, the real and imaginary parts are treated separately. The imaginary part is easier to compute once 
the integration over M 2 has been performed. We find 



7re 



d 3 P Q{-£ p )e{£ p+q ) 
(2tt) 3 



£p-\-q &p 



{£ p S(uJq — £p+q + £p) — £p+q${w q + Ep+q ~ £p)} 



(58) 



which is vanishing when oj q — is imposed. 

The real part of the r.h.s. of Eq. (^) is simplified by noticing that the exchange p^ p + q leads to the vanishing 
of the products of Heaviside step functions. Moreover, setting uj q — yields 



8 M 2Z iq 2 +d M 2U {2) S-2e 2 f 



d 3 P e(-£ p ) 



/' me 2 ki 



{1 + ffJ (2tt)3 (£ p+q - £ p ) (l + ff 2tt 2 



_q_ 

kp 



where 



1 _ l_ ( 

2 ~ 2x V ~ 4 



In 



1 + - 

1 T 2 



(59) 



(60) 



is the Lindhard function |lq| . The integration over M 2 from M — » oo to is easily computed. The expansion in 
powers of q finally gives 



Zi,m 2 {uo) 



2An 2 k F 



+ 0(M~ 2 ), U 2M * (u ) = ^-u 2 + 0(M~ 2 ) 



(61) 



The non- vanishing coefficient of Uq, rnkpe 2 /47T 2 , is the Thomas-Fermi mass giving the screening length in the electron 
gas. The other coefficient functions remain unchanged at the one-loop approximation. 



V. SUMMARY 



A new non-perturbative evolution scheme is developed in this paper for the electron gas with Coulomb interaction. 
No expansion in the electric charge is needed, the only approximation is the restriction of the evolution of the effective 
action to a given subspace of functionals, i.e. the consideration of a given set of 1PI functions with a given energy- 
momentum dependence. The evolution equations, a system of coupled equations for the 1PI functions are to be 
integrated out numerically from the initial condition given by the tree-level action functional. 

Since the evolution is performed by adjusting the amplitude of the quantum fluctuations there is no conflict between 
the successive elimination strategy of the renormalization group and the existence of a non-trivial Fermi surface. 

The unique feature of this scheme is that by truncating the effective action at higher order in the fields or in the 
derivatives one obtains systematically improved, non-perturbative approximations. Since the computation is carried 



11 



out in continuous space-time and the time can be chosen to be real this method offers an alternative for the usual 
numerical simulations of field theories on the lattice. 

The method provides the solution of the theory in the sense that physically important quantities are easy to express 
in terms of 1PI functions, the final result of the evolution equation. The only approximation, the truncation of the 
effective action in the power of the fields and their gradient has transparent interpretation and can be controled 
numerically. In fact, the maximal power of a given field in the effective action gives order of the 1PI functions 
retained in "turning on" the interactions. The order of the gradient expansion corresponds to the order of the Taylor 
expansion of the 1PI functions around the Fermi surface and p = lu = for fermions and bosons, respectively. 

We considered in this work the independent mode solution of the evolution equation only. This approximation 
reproduces the one-loop results of the usual perturbation expansion without relying Feynman graphs. In order to go 
beyond well known results one has to consider the system of coupled differential equations for the coefficient functions 
in the effective action and integrate them out numerically. This work, the numerical determination of the transport 
coefficients for the Coulomb gas is in progress. 
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APPENDIX A: FUNCTIONAL DERIVATIVES 



For the right hand side of the evolution equation we find the second functional derivative matrix of the effective 
action, r^ 2 ) [yb\ -0, u}. Since we expand the effective action in the fluctuations in order to identify the terms 0{d 2 ) we 
ultimately need the functional derivatives up to fourth order. 

The first four derivatives of the effective action ( |2l| ) written as V = J x (iplK x ip x + L x ), are 



r (D 
r (i) 



+Z Q dtu x dt5 x ^ y - Zi\7u x VS Xt y ( 



Y d t u x d t 6 x 



(Al) 
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where 



#W = iX { *\u x )d t - X[ n \u x )£ + U[ n \u x ) + ^Y {n \u x )(d t u x ) 2 + ^Y} n \u x ){Vu x ) 2 , 
4 n) = ^^ n) (ux)(ftux) 2 ~ \z { r\u x ){Vu x ) 2 - tf 2 (n) (u*) 

and the coefficient functions X, Y and Z are considered at the value u x in the integrands. 
The second derivatives, 



(A2) 



r ( 2 ) 
V>ii'/' B2 

p(2) 



r <2) 



r (2) 



Y d t u x d t 5 x , V2 + Y 1 'S7u x VS X} y 2 
Y d t u x dt8 x .y 2 + Y{Vu x VS Xiya 



+ 



u x,y 2 x,Vi 
fix, yi fix, V2 

J X 

(fj> x YV> + Zjp) d t u x (8 x,y 1 dtfi x ,y 2 + 8x,y 2 dtfix, Vl ) 

(i>li> x Yi 1 ' ) - Vu x {fix, Vl ^fix,v2 + fix, y2 Vtfix, yi ) 



fix, yi j 

i>x, 



+ / tylfaYo + z Q ) d t s x,y 1 dtfi x ,y 2 + / (4^Fi - Z x ) VS 

a: ,2/1 V(5 Xj y 2 , 



yield the non-diagonal, fluctuation-independent contribution to the self energy, 

E^" 
,„i = I 



(A3) 



(A4) 



with 



yu,ip _ yi>,u _ _ / t 4(1) 
Pi, Pa Pi.Pi ^— Pa— Pi P2 ' 

r"^' — v-iV t ,« _ 4(1) / 

Pl,P2 _ P2,Pl _ ^ 1 -pi-p 2 V y -Pl-P2> 



Pi ,P2 



B q ^ 



Pl ,P2 T—q—p l —p2 V 7 ?, 



(A5) 



and 



S Pi, P2 = A q 2) ~ y ou Pi oj P2 - Yipi -p2, 



(A6) 



and A being defined in Eq. (pq). The coefficient functions A, Y and Z are considered at uq in the gradient expansion 
formulae. 

In order to find the contributions of the fluctuations to the self energy up to the quadratic order we need the 
following functional derivatives, 



r 



(3) 
.(3) 

1p yi Uy^ liyg 



8 8 K^8 

<Jx,y 1 <Jx,y 3 I ^ x u x,y 2 
8x,y 2 8 x ,y 3 X, 



Y d t u x d t 8 XiVa + YiX7u x V8 x . 



{8x,y 3 dtS Xi y 2 + 8 X ^y 2 dt8 Xt y 3 ) 



8x,y 2 ! 



VU X 8 Xm (8 Xi y 3 V8 Xt y 2 + 8 X ,y 2 V8 X! y 3 ) 



Y d t 8 x ,y 3 dtS x ,y 2 +YiVS Xy y 3 V5 Xt y 2 
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r (3) 

U V2 U V3 



r (3) 

U Vl Uy 2 Uy 3 



r (4) 

^Vl ^V 2 u V 3 u Vi 



U yi Uy 2 Uy 3 tPy 4 



U V1 U V2 U V3 ^VA 



p(4) 

U V\ U V2 U V3 U VA 



m Sx,y 2 Sx,y 3 Ki 2) ip x + / 5 Xyyi Y^ 1) d t u x {S x , y3 d t 5 X! y 2 + d x ,y 2 d t S x , y3 ) tp x 
J x 

+ / 5 X! y 1 Y} 1 'VU X (5 X ,y 3 V5 X! y 2 +S X! y 2 'V5 X! y 3 )'lp X 
J X 



J x,yi 



YodtS x .y 3 d t S x .y 2 +Yi\7S x .y 3 X75 Xt y 2 
S x , Vl S x ,y 2 5 X! yMK^i; x + 4 3 )) + J (^V-^o (2) + 4 2) ) dtu 

X {^x,y 2 ^x,y^t^x,yi H~ $x ,y\$x ,y%dt$x ,y 2 ^x,y\ ^x,y 2 ^t^x.ys ) 

+ J x (^ x yI 2) -z[ 2) )vu x 

x (3x,y 2 dx,y 3 'V5 x ,y 1 + ^x,yi^x,y 3 V$x,y 2 + ^x,yi^x,y 2 ^^x,y 3 ) 

+ /(v4^o (1) + 4 1) ) 

J x 

x (S X:yi d t 5 X:y2 dt8 Xt y 3 + S xm d t S Xy y 3 dtS Xtyi + S x .y 3 dtS x . yi dtS x .y 2 ) 



(1) 



./a; ./a; 

^^a;,y4 H~ $x,y 4 ^^x,y^ ) 

/ ^x,yi^x,y 2 ^x,y 3 ^x,y 4 K^ ^1p x 
J x 

H~ / iftx^x, 3/4^0 dfU x ($x,y2^x,y3®t&x,yi ~l~ $x,y\$x,y%$t$x,y 2 ^x,yi $x,y 2 3t$x,y% ) 

H~ / Ipx^x, 3/4^1 {.^x,y2^x,y3^8x,yi "i" $x : yi$x, y^^^x, y 2 ^x,yi &x,y 2 ^^x,y% ) 

J a: 

x,yidtS x ,y 2 dtS x .y 3 + S x ,y 2 dt5 Xy y 3 dt5 Xyyi + 5 Xy y 3 dtS Xyyi dtS Xy y 2 ) 

J x 

+ / ip x 8 x , y4 Yi (S XyVl \76 Xy y 2 V6 Xy y 3 

~\~ dx,y 2 ^3 X ,y 3 ^S X ,yi ~\~ ^x,y 3 ^S Xy y 1 \7S Xy y 2 ) , 

</ X 

— / 3x,yidx,y 2 3x,y 3 ' i PxK x <5x,y 4 

/ ^x^x.yiY-^ Vu^ (<^x,i/2 ^x.93 V5 x ,9i &x,y\ &x,y 3 V5 X , y 2 + <5x,i/i ^x.92 ^^x,y 3 ) 

^ X 

— / Vx^x,j/ 4 ^o ' (^x, Vl dtS Xy y 2 dtS Xy y 3 + S Xyy2 dtS Xy y 3 dtS XyVl + 5 X y y 3 d t 5 X y y flfS X y y 2 ) 

J X 

/ V'x^x,94^ / i (^x,yi V(5 x .y 2 V(5a;,y 3 + ^x,92^^x,2/3 V5 x ,9i "I" &x,y 3 V$x,9i ^$x,y 2 ) ; 

/ 3x,y 1 8 x ,y 2 5x,y 3 3x,y4 (4^ 4 Vx + 4 4) ) 

J X 

+ / (V'xV'x^ + Z 3) ) 9tWxUx,j/2 5 x,i/3^x,94 5 ^x,9i 
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~~\~$x,yi $x,ys ^x^y^t^x.y^ ^x.y\^x,y2^x,y^t^x^y^ H~ $x,y\$x ,yo,^x ,yz®t$x ,y4 
~\~$x,yi $x,ys y^^^x, yi ~t~ ^#,2/1 ^,3/2 ^c,3/4^^e!3/3 ^£,2/2 ^£,2/3 ^^£,2/4 

X ,yiS X:Vi dtS Xj y 2 dtS Xt y 3 

~^~^x ,y2^x,V4&t&x ,y3^t^x .y\ ^x,y2^x,y3^t^x ,yi&t&x ~l~ ^x,y^x,ys^t^x,y2^t^x,y4 
~^^x,yi &x,y2 ^t^x ,ys&t3x ,y4 "t - ^x ,y^^x ,y4^t^x ,yi^t^x,y2^J 

+ j (i>llPxYx 2 ' > - Z[ 2 ^ \ S Xi y 1 S Xi y 4 VS X: y 2 VS xm +S X: y 2 S X: y i VS Xt y 3 VS X: y 1 

~^~^x,y2^x,y3^^x,yi^^x.y4 ~t~ ^a:, 3/1 ^£,2/3^^3?, 3/2^^,3/4 
~^~^x,yi ^x,y2 ^^x,y3 ^$x,y4 ~t~ ^x, 3/3 ^£,3/4^^2; ,3/i^^a:, 3/2 J ■ 

The matrix elements (pj) needed for the right hand side of the evolution equation can be read off easily, 



where 



'/— Pi— P2 



Pl,P2 



Pi ,P2 PI ,P2 



P2,Pl 



/ ^pi,p 2 ? 7-'?-pi-P2 r 7g> 



i/i+i 



Pl,P2 P2,Pl 



Pl,P2 P2,P1 



-9-P1-P2 



Pl,P2 



%tp! = / ^1,P2^^- 

91,92 

Fp 1 ,p2 1 lq*l J -q-p 1 -p2 — / Gp 1 , fp 2 rj qi ri q2 , il;_ qi _ q2 _p l _p 2 , 

J 91,12 

Hpi,P2V-pi-P2 + / JpllptVqi^i'- 

J 11 ,92 



-91-92-P1-P2 



(79 

Pl,P2 



I" / Jpl,p 2 2 ' q3 ' r lqi r lq2^Pq 3 ' l l J -qi-q2-q3-pi-P2 

'91,92,93 



^P2 } _ *W g (a;g + w Pl + lj P2 ) + Yig • (9 + pi + P2) 



^Pl,P2 7 ?9 7 ?-9-Pl-P2 7 



D q 

Pl,P2 



_ A (2) 



Pi ,P2 

F q 

Pl,P2 



r^9i 
Si 



9— Pl— P2 

(3) 



- Fo^g^p 2 - Fit? 'P2, 
., L " - ■/! -92-P1-P2 _ ^O^ [^91^92 + Hi + W 92H 2 ] - [ll ■ 92 + (qi + 92) -P2]j, 

Af} - >o^9 w P2 -YiQ-Ps, 

If = \ {4i } - Y ° 1] K + «>») + - i"/^ [P2 • (<Zl + 92 ) + 91 • 92]} , 

,P2 =^0 [( u pi + w p 2 ) 2 + w pi^p 2 ] + z i [(Pi+P2) 2 +Pi -Pa] +t r 2 3) > 



7-91,92 

Pl,P2 



791,92,93 
Pl,P2 



if 9 

Pl,P2 



= A 



(3) 

-9i— 92— Pi— P2 
I (4) 



- y o (1) Ki ( w pi + W P2 ) + w pi w p 2 ] - Y i X ' [91 • (pi + P2) + Pi ■ pa] , 



(i) 



2 I -91— 92 — 93-P1-P2 

,(2) 



F (2) [w 9l w g2 + + uj q2 )(uj pi + uj P2 ) + uj pi ajp 2 ] 



-*i [?1 ' 92 + (91 +92) • (pi +P2) +Pi -Pi] >, 



Zq 2) [{LO q + LO pi + L) p2 ) 2 + U q ((jJ pi + LU p2 ) + LU pi Up 



+ Z{ 2) [(q+ Pl +p 2 ) 2 -q-( P l+P2)-pi-P2] +/7 2 (4) |. 
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APPENDIX B: GRADIENT EXPANSION 



Here we provide more detailed expressions for the gradient expansion of the evolution equation Eq. (|17|). 

The fluctuations of the photon field to the evolution equation will be retaind up to the level O^ipr/ 2 ). These 
contributions from the expansion (^) are inserted on the right hand side of the evolution equation Eq. ([l?]). We 
introduce the notation 



■], B n = J2 T t \W<* 1 ,:.,a n ,J,j], 



(Bl) 



where the summation is over index sets where the number of value "fl" is 1 or 2. We find for the necessary terms 
-Tr[£W flj . Jt ] = 

Jp,p> 



c ^2 #' 



T 2 



Op<-r p \<T-p,p< <~r-p'L,-p< :P » Up" ^ 



-p",p 



0,p<D p ><J™ U p ,y,DpiiT, 



~P",P 



+Li -p,p' 1j p' l '-p',p" u -t >" a -p", P + a -p,p> u -p' a - P > , P " u p" ^-p" ,p + a -p,p' 1J p' cr -p',p" 1J p" lj -p",p 

^ a -p,p' L '~P' lj -p',p" 1J P" (J -p",p lj -p,p' 1J p" J -p',p" { ~'-p" a -p" ,p + ^ '-p,p> U P' -p> ,p" U P" ° -p" ,pJ ') 



C-p^pK V-ptf^-p' a -p> ,p" l*-p" ^-p» ip "> U P'" ^ 



-P"',P 



-(TZZ „/G- v <h, Y „, r .»-L>n"<J 



-p,p f p —p' 
tlpU 7-) ^uu 



p" w _pii pin Upin Zj 
t 



-ll/j-U 



+ 2j -P,P' ^ ^-P' ,P" U P' ,lj -p" ,p'» , 

for the electron contributions and 



p"',p ^-p,p l±J p " -p' , 
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for the photon induced terms. As discussed below Eq. (43) the sum of the terms listed in Eqs. (B2) and (B3) yields 
the evolution for Yq and Yi. 



1G 



